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ABSTRACT 



Using photometric catalogs from wide-field CCD images, we derive the cluster- 



o^ : 

Q\ , galaxy correlation amplitudes, Bgc, for 47 low-redshift Abell clusters. We apply a 

^~~^ \ number of tests to establish the robustness of the Bgc parameter as a quantitative 

^~ I measure of richness of galaxy clusters. These include using different galaxy luminosity 

p^H ' functions to normalize the excess galaxy counts, counting galaxies to different absolute 

Q^ ' magnitude limits, and counting galaxies to different cluster-centric radii. These tests 

^~~' . show that with a properly normalized luminosity function, the Bgc parameter is 

,-H . relatively insensitive (at better than 1/2 of an Abell Richness class) to magnitude 

^^ I limit, areal coverage, and photometric errors of up to about 0.25 mag. We compare the 

QQ I Bgc values to both the Abell Richness Class (ARC) and Abell count numbers {Na)- It 

^N ' is found that there is a good correlation between Bac and Na for Abell clusters with 

^^ ' z^O.l, with a dispersion of about one ARC; whereas the Abell Richness classifications 

5*^ ' for Abell clusters at z^O.l are much less well correlated with the true cluster richness. 

0\ • 

■^ ■ We also find evidence that the richness of ARC>3 clusters has the tendency of being 

f i . 

p . over estimated in the Abell catalog. 

^H ' To appear in The Astronomical Journal, May 1999 

Cd ■ 

> ■ 

"^ . 1. Introduction 

C^ ■ The richness is an important defining attribute of a galaxy cluster, yet has always proven 

to be difficult to quantify rigorously. Abell (1958), in creating the first comprehensive catalog of 
rich galaxy clusters, defined the richness classification which now bears his name based on the 
number of galaxies within 2 mag of the third ranking galaxy in an angular area of radius 1.7'/z, 
or about 3 Mpc (for Hq = 50 km s ~^ Mpc ~^). The procedure that Abell used in counting the 
number of galaxies from Palomar Sky Survey Prints must be considered approximate, and difficult 
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to extend to higher redshift. The accuracy of the Abeh classification for Abeh clusters is accurate 
probably to only about 1 class (e.g., see Dressier 1978). Uncertainties can arise from many 
sources: inaccurate estimate of photometry, poor star-galaxy classification, incorrect background 
correction, contamination from foreground or background clusters, and large uncertainty in the 
estimated redshift. In addition, the relatively large area used for the counting increases the chance 
of contaminating projections. At higher redshift, most of these problems would only be magnified, 
making a similar richness classification impractical without additional data such as accurate 
photometry and redshift. Furthermore, at redshifts higher than --^0.2, it becomes increasingly 
more difficult, if not impossible, to specify, for example, the third brightest cluster member, due 
to the large number of foreground galaxies. 

A robust richness measurement for galaxy clusters is an important observational quantity. 
Beside being one of the primary attributes of a cluster, it can be used as a baseline for comparing 
other properties of clusters in a meaningful way. This is especially important as advances in 
technology have allowed us to study galaxy clusters at different epochs and at increasingly higher 
redshifts, up to and above 1 (e.g., Dickinson 1997). Ultimately, a more fundamental parameter 
is the mass of the cluster. However, there is evidence that the richness of a cluster, properly 
determined, can be a good predictor of the mass of the cluster (Carlberg et al. 1996, and Yee et 
al. 1999) based on results from the CNOCl cluster redshift survey. 

Various attempts have been carried out to produce a more accurate measure of the richness 
of Abell clusters. These include the Nq^^ of Bahcall (1981), and the revised Abell Count Numbers 
,Ny{, of Abell, Corwin, & Olowin (1989, hereafter AGO). To arrive at a robust and rigorous 
quantitative measure of cluster richness, three observational issues must be taken into account 
when counting the number of galaxies in a cluster: the angular coverage, the depth of the 
observation, and the background correction. The availability of wide field digital imaging due to 
advances in large format CCDs has enabled us to address these issues quantitatively much better. 
The wide field allows us to sample the bulk of even very low redshift clusters, and digital imaging 
provides us with accurate photometry, which is vital for determining the depth of the galaxy 
counts, and correcting for background. In addition, a useful quantitative richness measure is one 
that can be applied uniformly to different samples at vastly different redshifts. 

The cluster-center galaxy correlation amplitude (Bgc), first used by Longair & Seldner (1979) 
to quantify the environment of radio galaxies, is one parameter that appears to be well-suited 
as a robust richness parameter. By assuming an universal luminosity function (LF) and spatial 
profile, the excess galaxy counts around any reference point can be properly corrected to create an 
unbiased estimator of the richness of the environment. Yee & Green (1987) extended the method 
to study the environment of quasars by using a "universal" LF with evolution that is self-consistent 
with the background galaxy counts. The ability to track the evolution of the LF at different epochs 
is extremely important in comparing richness of clusters at different redshifts in a meaningful way. 
However, although such a method has been used in many studies of the environments of active 
objects, fully digital photometric data have not been applied to the quantification and calibration 



of the richness of ordinary galaxy clusters, such as the Abell clusters, and to test the robustness of 
the parameter. 

In this paper, we present the richness measurements, in the form of Bgc, for a sample of 47 
Abell clusters based on digital photometric data. These measurements can be used as a calibration 
between the Bgc parameter and the more traditional Abell richness class (hereafter, ARC). We 
show that Bgc is a robust richness estimator, provided that a galaxy LF well matched to that 
found in clusters is used for normalization. We find that there is a good correlation between the 
Na and Bgc for clusters at z^O.l with a dispersion of about 1 richness class. However, z > 0.1, 
Na often over-estimates the true richness of Abell clusters. 

In §2 we present the sample and the data set. The method for deriving Bgc is briefly 
discussed in §3, along with the LFs used in the computation. In §4 we test the robustness of the 
Bgc parameter by investigating the dependence of the parameter on various quantities such as the 
LF, magnitude limit, and sampling area. Section 5 discusses the results and compare the ARC 
with the Bgc parameter. A summary is presented in §6. Throughout this paper, unless otherwise 
specified, we use Hq = 50 km s~^ Mpc ~^ and qo = 0.1. 



2. The Cluster Sample and Data 

The data used in this work are from a sample of Abell clusters originally chosen for a study 
of correlations between X-ray and optical properties (the Low-Redshift Cluster Optical Survey 
(LOCOS), see Lopez-Cruz &: Yee 1999). The detailed description of the sample, observation, 
and data reduction can be found in Lopez-Cruz (1997, hereafter LC97). Here, we present a brief 
summary. 

The clusters were selected from a compilation of bright X-ray clusters made by Jones & 
Foreman (1998) with the following criteria: (1) galactic latitude |6| > 30°; (2) 0.04 < z < 0.20; 
(3) ARC>1; and (4) the declination S > 24°. Although there is an imposed criterion of ARC>1, 
data for a number of ARC=0 clusters were also obtained, often for the purpose of filling Right 
Ascension gaps during the observing runs. The ARC=0 clusters are not considered as part of the 
fair sample. This is because they were not picked randomly; ones that appear to be rich or have 
well-defined X-ray emission were preferentially chosen. Nevertheless, they are analysed along with 
the other clusters, but are used only when such a selection effect does not affect the particular 
conclusion. As a comparison, the central region of Coma (A1656, z = 0.0232) was also observed. 
The total number of Abell clusters in the sample is 47, with 36 satisfying the criteria listed above. 
The cluster names, along with their redshifts and Abell Richness Classes, are listed in Table 1. 

Direct images in B, R, and I were obtained using the T2KA CCD camera at the 0.9m telescope 
at KPNO over five runs in 1992 and 1993. Photometry was calibrated to the Johnson-Kron-Cousin 
system of Landolt (1992). The CCD has 2048x2048 pixels with a scale of 0.68" pixel^i, giving a 
field of view of 23.2' x 23.2'. The integration times varied between 900 to 2500 second, depending 



on the filter used and the redshift of the cluster. 

Nine control fields at least 5° away from any cluster fields were also observed for the purpose 
of background count correction. These fields were observed and reduced in an identical manner as 
the cluster fields. 

Photometric catalogs of all objects in each field were created using the program PPP (Yee 
1991) to perform automatic object finding, photometry, and galaxy-star classification. The typical 
100% completeness magnitude in i? is -^ 21.5 mag. This provides a range of completeness in 
absolute R magnitudes at the redshift of the clusters from -15.5 to -18.5. 



3. Analysis 

3.1. The Bgc Parameter 

The cluster richness parameter Bgc is defined as the amplitude of the cluster center-galaxy 
correlation function: 



e(r) = BgcT-^. (1) 

This parameter was first used by Longair & Seldner (1979) to quantify the richness of radio 
galaxy environment. It was subsequently adopted by Yee & Green (1984) to study the environment 
of quasars using a self-consistent LF estimated from the same data. Further applications of the 
parameter to the study of the environments of quasars and radio galaxies can be found in Yee & 
Green (1987; hereafter YG87); Ellingson, Yee, & Green (1991); Prestage & Peacock (1988); and 
Yates, Miller, & Peacock (1989). This method has also been used to quantify the environment 
of BL Lacertae objects (Smith, O'Dea, & Baum 1995; Wurtz et al. 1997), Seyfert galaxies (De 
Robertis, Yee, & Hayhoe 1998), and the richness of a small sample of Abell clusters (Anderson & 
Owen 1994). Detailed descriptions of the method adopted for this study can be found in YG87, 
and also Andersen & Owen (1994), among others. Here, we provide a brief outline. 

We observe galaxies in projection on the sky, which allows us to measure the angular two-point 
correlation function of galaxies, uj{6), as a function of the angle on the sky. The function to (9) 
can be approximated by a power-law of the form (e.g., Davis & Peebles 1983): 

u;i9) = AggQ^-\ (2) 

where Agg is the galaxy-galaxy angular correlation amplitude. For our purpose, we can interpret 
a;(^) as the distribution of galaxies projected on the sky around a reference point, such as the 
center of a galaxy cluster. In this case, the amplitude, now notated as Ag^^ can be measured 
directly from an image by counting the excess (i.e., background corrected) galaxy counts, A„ei, 



up to a certain apparent magnitude around this reference point within some 9] i.e., if a fixed 7 is 
assumed, then Age = -jf^ ~^ 6"'"^ , where Ni,g is the background counts within 6. 

The Bgc amphtude can then be estimated via a deprojection of the angular correlation 
function into the spatial correlation function by assuming spherical symmetry, giving the relation 
between Bgc and Age (see Longair & Seldner 1979): 

where Ni,g is the background galaxy counts to apparent magnitude mo; ^[M{mQ^z)\ is the 
integrated LF of galaxies up to the absolute magnitude M corresponding to rriQ at the redshift 
of the cluster z. I^ is an integration constant arising from the deprojection, (with I^ = 3.78 for 
7 = 1.77), and D = -^{qoz + (go — l)[(2?o-2 + 1)2 — 1]/(7q(1 + z)^} is the angular diameter distance 

to z. 

As eqn (3) indicates, the computation of Bgc requires a knowledge of the luminosity function 
(which provides a normalization) and a background count correction. These are described in the 
next two subsections. 



3.2. The Control Field Counts 

The background counts used were obtained from the same cluster imaging program using the 
same telescope-CCD-filter combination and the same reduction procedure, which are described in 
detail in LC97. Briefly, 9 fields with integration times similar to those of the cluster fields were 
obtained at least 5° from any of the program clusters. The control field counts were combined 
using the technique described in Yee, Green, & Stockman (1986). Figure 1 shows the counts. 
These counts are in good agreement with the red band counts from other investigations (e.g., 
Driver et al. 1994, Metcalfe et al. 1991, and Yee et al. 1986). Using background counts obtained 
in an identical manner as the cluster data removes one of the major systematic uncertainties in 
the determination of Bgc- 



3.3. The Luminosity Function 

The choice of the LF used is a crucial one in the determination of Bgc- The sensitivity of 
the Bgc value to the LF will be discussed in detail in §4.1. The assumption of a universal LF is 
not strictly correct, as field and cluster galaxies have different population mixes and evolution 
histories. However, for the computation of Bgc using galaxies brighter than --^ — 18 mag, the 
varying faint-end slope of cluster galaxy LF (Lopez-Cruz et al. 1997) has little or no effect. We 
choose to test three LFs derived from very different data sources and methods. They are all 



represented functionally by the Schechter function (Schechter 1976). The three functions are (1) 
the LF derived from the clusters themselves (Lopez-Cruz & Yee 1999, hereafter the LCY LF); 
(2) the low-redshift field galaxy LF from King & Ellis (1985, hereafter the KE LF); and (3) the 
field galaxy LF from the CN0C2 redshift survey of intermediate redshift galaxies (Lin et al. 1999, 
hereafter the CN0C2 LF). The parameters of the LFs are listed in Table 2 and the LFs are 
described in the following. We will first discuss primarily the shape of the LFs used and defer the 
description of the normalization of the LFs to later in the section. 

We derive the shape of the LCY LF using the cluster data themselves. The derivation of 
the LF for this sample of clusters is described in detail in LC97, Lopez-Cruz et al. (1997), and 
Lopez-Cruz &: Yee (1999). The bright end is fitted to a Schechter function with a = —1.0. The 
result (LC97, Lopez-Cruz & Yee 1999) shows that the sample has < M* >= -22.26 ± 0.29, 
where the uncertainty is the dispersion of the distribution. However, cD clusters appear to have a 
larger range of M* which is apparently correlated with richness (LC97, Lopez-Cruz & Yee 1999). 
Removing the cD clusters produces a < M* >= —22.32 it 0.26. We adopt an "universal" LF for 
the clusters with M* = —22.3 (at a median redshift of 0.07) and a = —1.0. 

For deriving richness measurements for clusters at significant redshifts, the evolution of 
galaxies comes into play. Counting galaxies to a fixed absolute magnitude will in general over 
estimate the richness for high redshift clusters, as even mild passive evolution will mean effectively 
counting deeper into the LF. Hence, following YG87, we add a mild luminosity evolutionary term 
parametrized as —Qz mag, where Q ~ 1.4, to M* to compensate for this effect. (YG87 used a 
second order polynominal for M* evolution, we have simplified this to a single term in z, since 
any such measurement is currently extremely approximate.) There is good evidence that the LF 
in clusters evolves, to first order, in such a manner from galaxies associated with quasars (YG87), 
and from the CNOCl cluster redshift survey (Yee et al. 1999). Since the median redshift for the 
LC97 clusters is 0.07, we correct the M* to redshift to ~ 22.2. 

For a consistency check with previous applications of the B parameters in analyses of the 
environments of quasars (e.g., YG97), radio galaxies (Yates et al. 1989), BL Lacs (Wurtz et 
al. 1997), and Seyfert galaxies (De Robertis et al. 1998), we also compare Bgc values using the 
field galaxy KE LF with its separate components for different morphological types. The KE 
LF was derived as a preliminary result using a small number of low-redshift galaxies from the 
AAT/Durham redshift survey. Here, the M* values are transformed into the Cousin R band from 
the Gunn r band as listed in YG87. The same luminosity evolution as that for the LCY LF is 
incorporated into the KE LF. 

The final LF used is the interim LF derived from the CN0C2 field galaxy redshift survey (Yee 
et al. 1997, Lin et al. 1999). The LF is determined based on about 2300 galaxies with redshifts 
between 0.1 and 0.55, covering two patches on the sky. The galaxies are separated into early, 
intermediate, and late spectral types based on 5-color photometry. The LF for each component 
is then fitted to a Schechter function with two evolution parameters, P and Q, in the form of 



$(z,M) = $(0,M - Qz)W^-^^'' (see Lin et al. 1999 for details). Hence, the evolution of this LF 
is derived explicitly with Q measuring the luminosity evolution, and P, the density evolution (see 
Table 2). 

An important parameter in the LF for the determination of Bgc is the normalization constant 
(j)*. The value of Bgc is inversely proportional to it. As discussed in YG87, the best way to 
minimize systematic effects in Bgc is to determine a self-consistent LF with a normalization 
constant that reproduces the background counts as a function of apparent magnitude. We hence 
renormalized 0* by multiplying it by (j)' which is determined by fitting the modeled counts to the 
background counts. We determine (f)* for the LCY LF directly by minimizing the x^ between 
model counts generated by using the LF and the background count data. For the KE and CN0C2 
LFs, the relative 0*'s for each component are retained, but an overall normalization is applied 
based on the background count fitting. We found that the renormalization of the CN0C2 LF, 
which is derived using the same photometry band, is small, less than 10%. The modeled counts 
for the various LFs are shown in Figure 1. 



3.4. Results and Error Estimate 

The Bgc values obtained using the LCY LF are listed in Table 1. The values obtained using 
the other two LFs are similar, with variations of less than '^10%. A more detailed comparison of 
the dependence of Bgc on the LF is left to Section 4.1. For each cluster, we count galaxies down 
to a limit of M^ = —20 in a circle of 0.5 Mpc radius nominally using the brightest cluster galaxy 
(BCG) as the center of the cluster. However, we note that a few clusters have their BCG not 
sitting in the obvious center of the distribution of galaxies. The most significant deviation occurs 
in A168, where the BCG is almost isolated from the cluster. In this case, we have chosen the 
brightest galaxy in the high density clump (with the correct color) to be the center of the cluster. 
Choosing the BCG as the center for A168 would produce a Bgc almost 3 times smaller. 

The uncertainty for the Bgc parameter is computed using the formula: 

ABgc _ {Nnet + L32Ar,3)V2 



^gc ^^net 



(4) 



where N^et is the net counts of galaxies over the background of iVf,g. This is a conservatively large 
error estimate as it includes the expected counting statistics in Nnet and the expected dispersion 
in background counts. The factor 1.3^ is included to account approximately for the additional 
fluctuation from the clustered (and hence non-Poissonian) nature of the background counts (see 
LC97 and YG87). This is a fair estimate of the uncertainty of the richness of the cluster, as 
it includes the Poissonian uncertainty in the net counts, which is equivalent to the uncertainty 
in the number of galaxies drawn from a Schechter function with a fixed normalization (hence 
of fixed richness). However, we note that the inverse of this error estimate is an underestimate 



of the statistical significance of the excess galaxy counts in the sample, which would be better 
represented by Nnet/'^-^^i^^bg)^ ■ Hence, at these low redshifts, while the average uncertainty 
in the Bg^ parameter is about 25%, the significance of the excess counts over background is 
considerably better at about the 2Qa level. 



4. The Robustness of the Bgc Parameter 



In this section, we examine the results of exhaustive tests on the dependence of Bgc on various 
quantities, such as the LF shape and parameters, sampling and spatial limits, and photometric 
errors. These tests not only provide a clear indication of the robustness of the Bgc parameter, but 
also the pitfalls in applying the method. 



4.1. The Dependence of Bgc on the Luminosity Function 

A richness measurement based on the number of galaxies in a cluster is invariably tied to the 
assumption of the luminosity function of the galaxies. If clusters have greatly different galaxy 
LFs, then any such richness measurement would be rendered meaningless. In this section, we test 
the sensitivity of Bgc to the LF used in converting Age to Bgc- For these tests we vary the LF 
parameters but fix the sampling limit at M^ = — 20 mag. This sampling limit was chosen with the 
expectation that it minimizes the effect of having an incorrect LF. Varying the LF with different 
sampling limits causes more complex behaviors; these effects are tested in §4.2. 



4-. 1.1. Comparing Results from Different Luminosity Functions 

We have computed Bgc using three LFs derived using very different data and method. 
The LCY LF is based on the cluster LF measured from the same set of data, and describes 
cluster galaxies the best; while the KE and CN0C2 LF are based on field galaxies at 
different redshifts. The three give remarkably similar results with the mean of the ratio of 
Bgc{ljCY):Bgc{C^0C2):Bgc{KE) being 1.00:1.10:0.92, with the mean value for Bgc(LCY) being 
1246=b417, where the uncertainty is the dispersion of the Bgc values in the sample. This result 
demonstrates that the systematic uncertainty arising from the not-strictly-correct assumption of a 
"universal" LF for both cluster and field galaxies is of the order of 10%. The ratios indicate that 
the CN0C2 LF is slightly fainter overall compared to the LCY LF, while the KE LF is somewhat 
brighter. This indeed is the case when comparing the M*'s of early type galaxies in the two 
multi-component LFs with that of the LCY LF. 

For the remainder of the paper, we will use the Bgc values derived using the LCY LF as the 
fiducial set, and compare the results from the various tests to it. 



4-1-2. The Effect of an Incorrect Luminosity Function 

First, we test the effect of an "universal" LF with incorrect parameters. In varying the M* 
and a of the LF, we also adjust the normalization constant of the LF, so that it reproduces the 
control field galaxy counts up to the magnitude of sampling. Note that this is an important step 
to make sure that the Bg^s are obtained with a self-consistent LF model. This is different from 
testing the effects of statistical fluctuations of the LF of individual clusters from the mean for 
which one would not perform the renormalization (see §4.1.3). 

In these tests we first adjust the LCY LF by altering M* in steps of ±0.25 mag. Each time 
we recompute the normalization constant for the LF by fitting the background counts. We note 
that the fit to the background counts deteriorates with each step, but remains reasonable up to 
a AM of ±0.5 mag. The results are summarized in Table 3. The results for M* being off by 
±0.25 mag are quite acceptable, causing variations of about 10%, similar to the effects of using 
the three different LFs. This is not surprising, as the uncertainties in the three LFs are of the 
same order. Even when the error is 0.5 mag, the effect is about 20%, comparable to about 1/2 an 
Abell Richness Class (equivalent to ABgc of about 200; see Section 5). 

We next alter the faint-end slope a of the LF, again, with renormalization based on fitting 
the background. Here, we change a by steps of 0.15. The results are tabulated in Table 4. Again, 
one sees that changing the slope by as much as ±0.3 alters the Bgc values by less than 20%. 

These tests, along with the results from using three different LFs, demonstrate that a LF 
which self-consistently reproduces the background counts can be substantially incorrect and still 
produces results that are accurate to better than 1/2 an Abell class. 



4-1-3- The Effect of Fluctuations in the Cluster Galaxy Luminosity Function 

In this section, we test the effect of statistical fluctuations from the mean in the LF of the 
cluster galaxies. Here, we vary M* and a, but do not adjust the normalization constant of the LF. 
Note that in the case of M* , it has the same effect as a photometric error of the same magnitude 
in the cluster data. 

We use the same step sizes for M and a as in §4.1.2, and the results are listed in Tables 3 
and 4. Here, we see that the effects are significantly larger. Specifically, a fluctuation in M* (or 
equivalently, photometric error in the data) of about 0.4 mag will produce a 20% (1/2 Abell class) 
effect. However, a fluctuation of 0.2 mag would still be quite acceptable, producing uncertainties 
of only about 10%. Since cluster M*'s have a dispersion of about 0.2 mag, this indicates that the 
instrinsic dispersion of the LF from cluster to cluster does not affect Bgc signiflcantly. Similar 
fluctuations due to variations in a of up to 0.3 are also seen. 
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4.2. The Dependence of Bgc on the Limiting Magnitude 

The independence of Bgc values on galaxy counting limiting magnitude is an important issue, 
especially when different fields are observed to different depths. If the correct LF is chosen for 
the computation of Bgc and the LF is truly universal, then Bgc should subject only to counting 
statistical variations when different limiting magnitudes are used. A substantial systematic 
variation in Bgc with limiting magnitude can be taken as an indication that a LF with the 
incorrect shape has been chosen. Moreover, for the purpose of producing robust results, there is 
an optimal range of absolute magnitudes for galaxy counting. Counting to too bright a limit will 
produce a large uncertainty from counting statistics. Furthermore, one is at the mercy of small 
intrinsic variations of M* in individual clusters, or equivalently, small photometric zero point error 
in the galaxy photometry. 

On the other hand, a larger uncertainty may also result from counting too deep into the LF, 
as the background counts rise more rapidly than the flat part of the cluster LF. In addition, in 
counting beyond Mr ~ — 18 one may encounter large variations in the faint end slopes of the 
cluster LF (Lopez-Cruz et al. 1997 and LCY), producing additional fluctuations. Ideally, one 
would like the richness parameter of a galaxy cluster to be representative of the total mass of the 
cluster, which, under the light-traces-mass scenario (e.g., Tyson & Fischer 1996; Carlberg, Yee, 
& Ellingson 1997), should be proportional to the total light. Hence, the richness of a cluster is 
best represented by the number of iy*-like galaxies, rather than by the number of dwarf galaxies; 
consequently, a robust measure of Bgc is best obtained using a magnitude limit corresponding to 
the flat part of the LF: from about 1 mag to 3 mag past M* . 

The data in our sample of Abell clusters sample to at least Mr = —18.5 mag for all the 
clusters; hence, they provide a good test on the robustness of Bgc versus counting magnitude, and 
the effects of varying the LF. For these tests, we compute the Bgc parameters using Mr = —22.0, 
-21.0, -20.0, -19.0, and -18.0 for comparison with the standard limit of -20 (i.e., about 2 mag past 
M*). The five derivations provide ratios of average Bgc of: 1.04:1.05:1.00:1.00:1.07. A comparison 
of the individual BgcS using sampling limits of -21.0 and -20.00 is shown in Figure 2. This test 
shows very strongly that, as expected, when a LF closely resembling that of the cluster galaxy LF 
is used, the sampling magnitude limit has little or no systematic effect on the derived richness. 
The variations in individual Bgc values are the order of the computed uncertainties, and are due 
to counting statistics. We note that the largest deviation occurs with the deepest sampling limit. 
This is expected as the steepening of the galaxy LF for some clusters begins at about Mr = —18, 
which would produce a larger count relative to the flat (a = 1) LF used. 

Next, we consider the effects of the combination of LF variations and more extreme sampling 
limits. We focus on testing the effect of incorrect LF and photometric errors when sampling to only 
about M*, a situation that may arise often when data which are not quite deep enough are pressed 
into use. The results are shown in Table 3. In general, at a sampling limit of Mr = —22 mag 
(about M* ) , the errors committed when M* is varied with renormalization of the LF (equivalent 
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to using an incorrect LF) are about twice as large as those when the samphng is to ~2 mag past 
M* (about -20 mag). However, when no renormahzation is performed, which is also equivalent to 
photometric uncertainty, the errors incurred when counting to -22 mag can be as large as 4 times 
those to -20 mag. For example, a +0.5 mag change in M* produces Bgc values 2.2 times larger 
when counting to -22 mag, compared to 1.3 times larger when counting is done to -20 mag. We 
note that the results for changing a are not nearly as drastic, as that affects mainly the faint end 
of theLF (see Table 4). 

We hence conclude that counting to about M* produces results that are extremely sensitive 
to errors in the assumed LF, photometry, and in general statistical fluctuation, with results that 
could be off by a factor of 1.5 to over 2 when the error in magnitude is 0.25 to 0.5 mag. This may 
explain the large discrepancies noted in the Bgc values of BL Lac objects between Smith, O'Dea, 
& Baum (1995) and Wurtz et al. (1997), as noted by Wurtz et al.. We also note that this may also 
be the possible explanation of the discrepancies of between the richness derived for Abell clusters 
by Andersen & Owen (1994) and this work (see §5.1). 



4.3. The Dependence of Bgc on Sampling Area 

Ideally, a robust richness parameter is one such that the counting area used should not have 
a significant effect on the richness derived. Such a parameter would allow comparison of the 
richness parameter from data taken with a wide variety of areal coverage. This is strictly possible 
only if the spatial distribution function of cluster galaxies has a similar form for all clusters. If 
an incorrect slope (7) for the correlation function is used, the Bgc value will be dependent on 
the sampling area. In the computation of Bgc it is assumed that the galaxy distribution drops 
off in the same fashion as the spatial correlation function for field galaxies. This assumption is 
adopted so that comparisons over a wide range of richness and environments, including that of 
field galaxies, can be made. However, this assumption is most likely incorrect for rich clusters, as 
various investigations have shown that the power-law slope of their galaxy density distribution is 
probably greater than 2 (e.g., Lilje & Efstathiou 1988, Pebbles 1993), somewhat steeper than that 
of the clustering of field galaxies. 

The effect of using an inexact 7 can be examined by comparing the Bgc values for our sample 
using different limiting radii. We have computed Bgc values (with the LCY LF) using sampling 
radii of 0.25, 0.5, 1.0 Mpc, all with a Mji sampling limit of -20 mag. The ratio of the mean BgcS 
for the 3 radii are: 1.03:1.00:0.94. The results for individual clusters with sampling radii of 0.5 
and 1.0 Mpc are shown in Figure 3. 

The results indicate that a change of a factor of 2 in the sampling radius causes a scattering 
of about 5% in the richness parameter. This small effect is reassuring, and provides confidence 
when comparing data sampled over metric field sizes of a factor of several, as in those used by 
De Robertis et al. (1998) for nearby Seyfert galaxies. The decreasing trend in the Bgc with larger 



12 



radii is consistent with existing evidence that cluster galaxy distribution is likely to be steeper 
than the slope of the field galaxy correlation function. However, since the effect is very small, 
we have chosen to retain the canonical 7 = 1.8 so that direct comparisons with field and group 
clustering amplitudes can be made. 



5. Discussion 

The various tests in Section 4 have established the robustness of the Bgc parameter as a 
measure of cluster richness. In this section, we compare the Bgc parameter with the Abell Richness 
Class, establishing a calibration between the two. We also discuss the correlation between Bgc and 
the Abell count number, which allows us to examine the accuracy of the Abell count number as 
a measure of richness. Based on these comparisons, we suggest a revised richness classification 
scheme based on the calibration between the Bgc parameter and ARC. Finally, we illustrate 
the advantage of using a robust richness parameter in the investigation of cluster properties by 
examining the correlation between richness and velocity dispersion. 



5.1. Comparison with Abell Richness Class 

One of the main goal of this investigation is to establish a calibration between a more 
quantitative richness parameter, such as Bgc, and the traditional Abell Richness Class. In Figure 
4, we plot the Bgc parameters versus ARC. We see in general a very broad correlation between 
Bgc and ARC. We note that because we have specifically chosen to image ARC clusters which 
appear to be the richest, they should not be considered in the correlation. 

Andersen &; Owen (1994) produced a qualitatively very similar plot of the two quantities 
using a different sample and very different type of data. However, quantitatively, there is a 
significant discrepancy in the scaling of the two quantities: Their Bgc values are close to a factor 
of two lower than ours (after correcting for Hq), with their average Bgc for ARC 1 clusters being 
~ 600 Mpc^^'^, compared to about 1000 Mpc~^'^ in our work. This discrepancy is perhaps not 
surprising, given that: (1) Andersen & Owen use POSS photographic data which may not be well 
calibrated, (2) their object counting goes down to only ~ M* , and (3) they did not use a LF that 
has been renormalized to the background. All these can contribute substantially to the systematic 
uncertainty, up to a factor of 2, as demonstrated by our tests in §4. 

Since the ARC is an approximate and quantized measure, we will next look at the comparison 
with actual Abell counts (Nj^) before revisiting the correlation between Bgc and ARC in §5.3. 
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5.2. Comparison with Abell Count Numbers 

A more detailed comparison between the ARC and the Bgc measurement can be obtained by 
using the Abell number counts (Na) as published by AGO. We plot Na vs Bgc in Figure 5. Taken 
as a whole, again we see a definite correlation between Bgc and Na but with a large dispersion. 
Specifically, it is noted that at large Na the correlation breaks down in that their corresponding 
Bgc values are smaller than indicated by Na- We also note that all the clusters in our sample 
with large Nas are those at higher redshifts. We hence break the sample into two groups, dividing 
them at z = 0.09. To increase the z > 0.09 sample, we also add 3 clusters from the CNOCl survey 
(Yee et al. 1999) which are AbeU clusters at z ~ 0.2. These are A2390, MS0451+02 (A520), and 
MS0906+11 (A750). 

Removing the high redshift sample provides a much improved correlation bewteen Bgc and 
Na- For the lower z sample, though the scatter is about one Abell class, there are few or no 
outliers. For the higher z sample, there is almost no correlation between Bgc and Na- The 
discrepancies come entirely in the high Na clusters which all have too large a Na for their B ge- 
lt is interesting to note that A665, which is the only Abell cluster classified as ARC 5, with 
Na = 321, has a Bgc value which indicates that its richness is overestimated by a factor of 3 by 
Na- a number of Abell clusters in fact have similar richness to A665, including A2390 and A401. 
Hence, we conclude that at z^O.l, the AGO counts provide reasonable estimates of the richness of 
the clusters; however, at zk.0.1, and/or for Na^^^O, the Na counts almost always over predict the 
cluster's richness. This is not a surprising result, as many authors have concluded that the Abell 
catalog becomes incomplete at zk,0.1 (e.g., Southerland 1988; Struble & Rood 1991b). 



5.3. A Bgc Calibrated Richness Classification 

Using Figures 4 or 5, we can relate the more robust Bgc measurements directly to the Abell 
Richness class, in essence, calibrating the classical richness class designations. The mediam 
Bgc value of ARC 1 clusters in our sample is very close to 1000 Mpc^'^. Hence, we propose to use 
1000 Mpc^'® as the fiducial Bgc value for the average ARC 1 counts. 

Abell (1958) and ACO used a non-linear relationship between the Abell counts and the 
richness class with the higher ARC having disproportionately higher counts, presumably to 
account for the few apparently very high Na clusters . However, §5.2 demonstrates that the 
very high Na clusters are probably much less rich. If we scale the Bgc values corresponding to 
each ARC counts using the fiducial of 1000 Mpc^'® as ARC 1, we see that there is a significant 
mismatch at ARC>3, as illustrated by the dotted line in Figure 4, with ARC 5 being equivalent 
to Bgc^ 7000. Yee et al. (1999) also found that none of the 16 clusters from the CNOCl survey, 
which contains the most X-ray luminous clusters from the EMSS survey (see Gioia & Luppino 
1994), have Bgc larger than --^ 2500 Mpc ~^'^. Hence, there appears to be no known clusters that 
approach the original definition of ARC 5. In view of this, we propose a linear richness scale 
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where adjacent richness classes change by the same increment of 400 Mpc^'^ in Bgc- This scale 
would cover most of the range of known cluster richness (including very rich, X-ray luminous, 
higher redshift clusters, such as those in the CNOCl sample), with the first 3 classes (0, 1, and 
2) matching the original ARC definition reasonably well. This revised classification definition is 
shown graphically in Figure 4. For Hq = 50 km s^^ Mpc ^^, the calibration for ARC to 5 is 
then: 600±200, 1000±200, 1400±200, 1800±200, 2200±200, and > 2400 Mpc^-^, respectively. 

We note that this calibration is substantially different from those provided by Prestage & 
Peacock (1988) based on a very small number of 4C radio galaxies situated in Abell clusters. Their 
cahbration of ARC and 1 having Bgc of ~380 and 850 Mpc^'^ is about 1/3 to 1/2 of an Abell 
class too high. Hence, previous descriptions of the Abell class equivalence for the environments of 
radio galaxies and quasars (e.g., YG87) need to be revised correspondingly downward. 



5.4. The Correlation Between Velocity Dispersion and Richness 

So far we have shown that a richness parameter such as Bgc, which uses photometric data and 
normalized by the LF and spatial distribution of galaxies, produces a robust measurement of the 
richness of a cluster. One would then expect that such a robust parameter will provide a better 
description of the clusters, and hence produce more meaningful correlations with other cluster 
properties. 

We will show, as a demonstration of the improvement expected, an example of such a 
correlation. One fundamental cluster property is the mass of a cluster, which can be estimated 
by the velocity dispersion (a). In Figure 6 we plot the correlation of a versus Bgc- The velocity 
dispersion measurements are taken from the compilation of Struble & Rood (1991a), and Fadda et 
al. (1996). We include only those clusters with 10 or more velocity measurements. There is a clear 
correlation between a and Bgc, and the correlation is much tighter than using A^^ as the variable 
for richness which is shown in Figure 7. Girardi et al. (1993) showed a similar improvement in the 
correlation between velocity dispersion and richness when Bahcall's (1981) A''o.5 parameter is used 
instead of the Abell counts. This improved correlation lends support that Bgc is a more robust 
richness measurement. 

The (T-Bgc correlation derived from the CNOCl cluster redshift survey (Yee et al. 1999) 
provides an even tighter correlation between the two quantities. This is almost certainly due to 
the fact that the velocity dispersion measurements from CNOCl are derived from a homogeneous 
set of velocity data with a large number of velocities for each cluster. We note that the 
(T-Bgc relationships derived from the two different samples (this work and CNOCl) are identical 
within the uncertainty. The cr-Bgc relationship from CNOCl is plotted as a dotted line in Figure 
6. Given the different photometric systems used (Kron-Cousin R vs Gunn r), the different redshift 
range of the sample, and the different instrumentations used for these two samples, this again 
points to the robustness of the Bgc measurements as a richness estimator. 
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6. Summary 

We have used a sample of 47 Abell clusters with photometric catalogs from CCD images 
to test the robustness of Bgc as a parameter for measuring the richness of galaxy clusters. We 
first tested the effect of the choice of galaxy LF, which is used to normalize the excess galaxy 
counts. Using three LFs derived from different sources, including both cluster and field galaxies, 
we found that as long as we adjust the density normalization constant of the LF so that it 
reproduces the background counts used to derive the excess counts, different choices of the LF 
produce variations in Bgc of approximately 10%. Furthermore, the parameters for the LF can be 
substantially incorrect (e.g., with M* being off by 0.5 mag), errors in Bgc can also be minimized 
by renormalizing the LF to fit the background counts. 

We also found that if a correct LF is used, counting to different absolute magnitude limits 
produces systematic variations of only a few percent. However, to minimize uncertainties 
introduced by LF discrepancy and photometry errors, the results are most reliable when the 
counting is done to about 1 to 2 mag past M* . 

Similarly, the canonical 7 = 1.8 for the spatial correlation function of field galaxies provides 
reasonably stable Bgc values as one varies the radius of the counting aperture. Altering the 
counting radius from 0.25 to 1.5 Mpc produces systematic changes in Bgc of less than 10%. 

We have compared our Bgc richness parameters with the Abell Richness Class and the Abell 
Count Number {Na_)- Limiting the sample to either ARC<2 or z^O.l, we found a good correlation 
between A^^ and Bgc- However, clusters designated as ARC>3 appear often to have their richness 
over estimated, sometimes by as much as a factor of 3. Since most of the ARC>3 clusters are at 
z > 0.1, this suggests that ARC classifications for clusters with z^O.l is not a good indicator of 
their true richness. We suggest a revised richness classification which is based on a linear scale in 
Bgc- Finally, we demonstrate the improvement in the correlation between richness and velocity 
dispersion when a robust richness parameter, like Bgc, is used. 

We wish to thank KPNO for the generous allotment of telescope time and assistance in 
carrying out the Low-redshift Cluster Optical Survey (LOCOS). HY acknowledges NSERC 
for financial support in the form of an operating grant. OLC was supported in part by an 
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from the Department of Astronomy and the School of Graudate Studies of the University of 
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Fig. 1. — Background galaxy counts in Cousin R from 9 KPNO O.Om fields. The dashed, dotted, 
and solid lines represent the count models created using the LCY LF, KE LF, and CN0C2 LF, 
respectively. 

Fig. 2. — Comparison of Bgc values computed using counting magnitude limits of -20.0 and -21.0. 
The former cut-off is equivalent to M* — 1.7. The dotted line is not a fit, but a line of slope unity. 

Fig. 3. — Comparison of Bgc values computed using counting radii of 0.5 Mpc and 1.0 Mpc. The 
dotted line is not a fit, but a line of slope unity. 

Fig. 4. — The correlation between Abell Richness Class and the Bg^ parameter. The median 
Bgc value for each ARC is shown as an asterix. The dot-dashed line indicates the expected 
Bgc values for each ARC which describes the relative median richness of ARC 1 and 2 well, but 
fails significantly for richer classes. The six bars connected by the dashed line illustrates a linear 
richness classification which matches the real richness of the various ARC reasonably well. 

Fig. 5. — The correlation between Abell count number and the Bgc parameter. Solid symbols 
represent cluster with z < 0.09; while open symbols mark clusters with z > 0.09. The solid star 
symbols are Abell clusters at ~ 0.2 from the CNOCl survey. The dashed line shows the best fitting 
line for the z < 0.09 subsample force-fitted through the origin, while the dotted line is the best 
linear fit for the z > 0.09 subsample. 

Fig. 6. — The correlation between velocity dispersion and cluster richness as parametrized by Bgc- 
The dashed line is the best linear fit; while the dotted line shows the relationship derived from 15 
CNOCl clusters (Yee et al. 1999). 

Fig. 7. — The much poorer correlation between velocity dispersion and cluster richness as 
parametrized by the Abell count number (Na)- 
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Table 1. Abell Clusters 



Cluster 


y 


ARC 


Bgc 


± 


A21 


0.0946 


1 


1509 


241 


A84 


0.1030 


1 


1000 


204 


A85 


0.0518 


1 


757 


174 


A98 


0.1043 


3 


1921 


268 


A154 


0.0638 


1 


1393 


230 


A168 


0.0452 


2 


973 


192 


A399 


0.0715 


1 


1449 


235 


A401 


0.0748 


2 


2242 


286 


A407 


0.0472 





1368 


225 


A415 


0.0788 


1 


534 


159 


A514 


0.0731 


1 


946 


196 


A629 


0.1380 


1 


1205 


224 


A646 


0.1303 





858 


196 


A665 


0.1816 


5 


2186 


290 


A671 


0.0491 





1266 


217 


A690 


0.0788 


1 


701 


175 


A957 


0.0437 


1 


1075 


201 


A1213 


0.0469 


1 


931 


189 


A1291 


0.0530 


1 


1122 


207 


A1413 


0.1427 


3 


1686 


257 


A1569 


0.0784 





768 


181 


A1650 


0.0845 


2 


1861 


263 


A1656 


0.0232 


2 


1242 


282 


A1775 


0.0700 


2 


1018 


202 


A1795 


0.0621 


2 


1430 


233 


A1913 


0.0530 


1 


954 


192 


A1983 


0.0430 


1 


974 


192 


A2029 


0.0768 


2 


1736 


255 


A2244 


0.0997 


2 


1698 


254 


A2255 


0.0800 


2 


2296 


289 


A2256 


0.0601 


2 


2174 


281 


A2271 


0.0568 





645 


164 


A2328 


0.1470 


2 


1941 


273 


A2356 


0.1161 


2 


944 


201 


A2384 


0.0943 


1 


1575 


245 


A2399 


0.0587 


1 


906 


190 



20 



Table 1 — Continued 



Cluster 


z 


ARC 


Bgc 


± 


A2410 


0.0806 


1 


699 


175 


A2415 


0.0597 





902 


190 


A2420 


0.0838 


2 


1230 


220 


A2440 


0.0904 





1053 


207 


A2554 


0.1108 


3 


1218 


222 


A2556 


0.0865 





828 


187 


A2593 


0.0421 





1211 


212 


A2597 


0.0825 





665 


172 


A2626 


0.0573 





945 


193 


A2657 


0.0414 


1 


740 


170 


A2670 


0.0761 


3 


1771 


257 



Table 2. Luminosity Functions 



LF 


galaxy type 


<!>' 


r 


Q 


M* 


Q 


P 


LCY 


cluster 


0.00226 


1.0 


-1.00 


-22.20 


1.4 


0.0 


KE 


E/SO 


0.000346 


1.79 


-1.00 


-22.45 


1.4 


0.0 




Sab 


0.000346 


1.87 


-1.00 


-22.14 


1.4 


0.0 




Sbc 


0.000346 


1.87 


-1.00 


-21.81 


1.4 


0.0 




Scd/Sdm 


0.000346 


3.24 


-1.00 


-21.29 


0.7 


0.0 


CN0C2 


early 


0.92 


0.00212 


0.09 


-21.52 


1.94 


-2.02 




intermediate 


0.92 


0.00111 


-0.69 


-21.84 


0.62 


0.59 




late 


0.92 


0.00053 


-1.45 


-21.44 


0.88 


2.68 
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Table 3. Fractional Change to mean Bgc Adjusting M 



Sampling Limit LF Renormalization AM*: -0.50 -0.25 0.00 +0.25 +0.50 

-20 yes 

-20 no 

-22 yes 

-22 no 



1.24 


1.11 


1.00 


0.91 


0.84 


0.80 


0.89 


1.00 


1.14 


1.32 


0.90 


0.95 


1.03 


1.18 


1.42 


0.58 


0.76 


1.03 


1.47 


2.24 



Table 4. Fractional Change to mean Bgc Adjusting a 



Sampling Limit LF Renormalization Aa: -0.30 -0.15 0.00 +0.15 +0.30 

-20 yes 

-20 no 

-22 yes 

-22 no 



0.92 


0.96 


1.00 


1.10 


1.19 


0.73 


0.86 


1.00 


1.14 


1.29 


1.39 


1.19 


1.03 


0.95 


0.87 


1.11 


1.07 


1.03 


0.99 


0.94 
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